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I.  INTRODUCTION 


A,  MEDITERRANEAN  SEA 

The  Latin  word  “Medius  Terrae,”  which  means  “inland”  or  “in  the  middle  of 
land,”  is  known  as  the  origin  of  the  name  of  the  Mediterranean  Sea.  The  Mediterranean 
region  is  surrounded  by  the  Europe,  Asia  and  Africa  continents  (Tanabe,  2013).  The 
Mediterranean  region  countries  are  seen  in  Figure  1 . 


Figure  1 .  Mediterranean  Sea  and  Surrounding  Countries 
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Adapted  from:  WorldAtlas,  2015:  Map  of  the  Mediterranean  Sea.  Accessed  24  August 
2015.  [Available  online  at  http://www.worldatlas.com/aatlas/infopage/medsea.htm] 


The  Mediterranean  Sea  stretches  over  an  area  of  2.5  million  square  kilometers, 
derived  from  nearly  3000  kilometers  in  longitude  and  1500  kilometers  in  latitude.  It  is  a 
partly-encircled  basin  that  is  linked  with  the  Atlantic  Ocean  through  the  Strait  of 
Gibraltar;  linked  with  the  Black  Sea  through  the  Dardanelles  Straits  and  the  Turkish 
Bosphorus.  The  Mediterranean  Sea  is  divided  by  the  Sicily  Channel  into  two  parts.  They 
are:  Eastern  Mediterranean  and  Western  Mediterranean  (Soto-Navarro  and  Criado- 
Aldeanueva,  2012). 
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The  Mediterranean  Sea  extends  between  the  30°  and  46°  north  latitudes;  5°50' 
west  and  36°  east  longitudes  (Boxer,  2014).  The  Western  Mediterranean  includes  the 
Alboran,  Balearic,  Tyrrhenian  and  Ligurian  Seas.  The  Eastern  Mediterranean  includes  the 
Levantine,  Adriatic,  Ionian  and  the  Aegean  Seas  (Theocharis  et  ah,  1998).  Major  parts  of 
the  Mediterranean  Sea  are  seen  in  Figure  2. 

Figure  2.  Major  Parts  and  Major  Sub-basins  of  the  Mediterranean  Sea 


Source:  Wikipedia  contributors,  2015:  Mediterranean  Sea.  Accessed  25  August 
2015.  [Available  online  at  https://en.wikipedia.org/w/index.php?title=Mediterranean 
_Sea&oldid=695 192772] 
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The  Mediterranean  Sea  is  a  quasi-eneireled  sea  nearly  landlocked  by  continents. 
When  correlated  to  the  other  oceans’  circulation  and  residence  time,  that  is  about  200- 
1000  years,  the  Mediterranean  Sea  is  known  by  its  short-lived  circulation  and  residence 
time,  that  is  about  70  years  (Group  et  ah,  2011). 

The  overall  length  of  the  Mediterranean  Sea  basin  is  3,860  kilometers,  and  the 
utmost  width  is  1,600  kilometers.  The  whole  area  enclosed  by  the  basin  is  nearly 
2,536,000  kilometers  square  and  the  entire  volume  is  almost  3,750,000  kilometers  cubic. 
The  Mediterranean  Sea  has  only  one  link  to  the  Atlantic  Ocean  through  Strait  of  Gibraltar 
channel  that  is  14  kilometers  wide  and  600  meters  mean  depth  (Papadopoulou  et  ah, 
2011).  Topographic  and  bathymetric  map  of  the  Mediterranean  Sea  is  seen  in  Figure  3. 


Figure  3.  Topographic  and  Bathymetric  Map  of  the  Mediterranean 

Sea 


Source:  Universita  degli  Studi  di  Pavia,  2005:  The  Mediterranean  Sea.  Accessed  02 
September  2015.  [Available  online  at  http://www-3.unipv.it/cibra/edu_Mediterraneo_uk] 
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The  average  depth  of  the  Mediterranean  Sea  is  nearly  1450  meters;  henee  it  is  a 
shallow  basin  when  it  is  compared  to  the  mean  depth,  almost  3850  meters,  of  all  oceans 
in  the  world.  While  the  deepest  depth  is  nearly  3733  meters  in  the  Western  Mediterranean 
Sea,  the  estimated  mean  depth  of  it  is  almost  1612  meters.  The  deepest  depth  in  the 
Eastern  Mediterranean  basin  is  5121  meters  (Papadopoulou  et  ah,  2011).  Depth  zones  of 
the  Mediterranean  Sea  are  seen  in  Table  1. 


Table  1 .  Depth  Zones  by  Area  of  Coverage  of  the  Mediterranean  Sea 


Depth  Zone  (meters) 

Area  (kilometers  square) 

0-200 

578000 

200-1000 

720000 

1000-2000 

230000 

2000-3500 

951000 

Over  3500 

57000 

Total 

2536000 

Adapted  from:  Papadopoulou,  K.  N.,  V.  Markantonatou,  and  C.  J.  Smith,  2011:  Options 
for  Delivering  Ecosystem-based  Marine  Management.  The  Mediterranean  Sea: 
Additional  information  on  status  of  threatened  ecological  characteristics  relevant  to  the 
Marine  Strategy  Framework  Directive,  30  pp,  [Available  online  at 
https://www.liverpool.ac.uk/media/livacuk/odemm/images/ 
MediterraneanSeaBackground.pdfj 


B,  OBJECTIVE 

Improvements  in  the  political,  military  and  economic  fields  in  the  Mediterranean 
region,  specifically  in  the  Eastern  Mediterranean,  requires  a  deeper  understanding  of  the 
hydrodynamics  structure  of  the  Eastern  Mediterranean  Sea,  along  with  its  natural  aspects. 
Thus,  the  spatial  and  temporal  variations  in  the  Eastern  Mediterranean  Sea  occupy  a 
critical  importance,  notably  in  terms  of  the  military  based  operational  objectives.  This 
thesis  will  scope  on  analyzing  the  data  obtained  from  World  Ocean  Database.  Analysis 
will  include  inter-annual,  seasonal  and  monthly  salinity,  temperature  and  current 
variation  fields  from  1960  to  2013.  These  analyzed  salinity,  temperature  and  current 
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profiles  will  lead  to  determine  the  general  current  patterns  and  thermohaline  structures. 
After  analyzing  the  data,  results  will  be  compared  with  the  previous  studies.  Thus,  a 
further  understanding  of  the  general  circulation  patterns  of  the  Eastern  Mediterranean  Sea 
will  be  presented. 

C.  THESIS  ORGANIZATION 

In  Chapter  II,  some  crucial  aspects  of  the  Eastern  Mediterranean  and  Eevantine 
Sea  hydrodynamics  will  be  explained.  They  are;  water  masses,  wind  fields,  nutrient  fields 
and  oceanographic  features  of  the  Eastern  Mediterranean  and  Eevantine  Sea.  In  Chapter 
III,  information  about  the  dataset  used  in  the  thesis  will  be  touched  on.  The  source  of  the 
dataset,  methods  and  instruments  used  for  collecting  the  data  will  be  mentioned.  In 
Chapter  IV,  methodology  for  creating  the  fundamental  principles  of  the  study  will  be 
explained.  In  Chapter  V,  analyses  of  the  data  will  be  conducted.  By  conducting  the 
analyses,  it  will  be  divided  into  sub  sections  for  deeper  understanding.  It  will  mainly 
focus  on  the  total  mean  and  seasonal  variations  for  the  temperature  and  salinity  profiles 
as  well  as  the  current  profiles.  In  Chapter  VI,  inter-annual  variability  of  the  Eastern 
Mediterranean  Sea  and  first  three  EOE  modes  will  be  studied.  In  Chapter  VII,  results  that 
are  taken  by  this  study  will  be  summarized  and  any  recommendations  for  further  studies 
in  the  future  will  be  pointed  out. 


5 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


6 


II.  OCEANOGRAPHIC  FEATURES  OF  THE  EASTERN 
MEDITERRANEAN  AND  LEVANTINE  SEA 


A,  EASTERN  MEDITERRANEAN  AND  LEVANTINE  SEA 

The  Eastern  Mediterranean  Sea  ineludes  four  important  major  basins  and  seas. 
They  are;  Ionian  basin,  Levantine  basin,  Aegean  Sea  and  Adriatic  Sea  (Alhammoud  et 
ah,  2005).  Major  sub-basins  of  the  Eastern  Mediterranean  Sea  are  seen  in  Eigure  4.  The 
Ionian  basin  on  the  west,  north-eastern  Africa  on  the  south.  Middle  East  on  the  east, 
Cretan  Archipelago  and  Asia  Minor  on  the  north  surround  the  Levantine  basin 
(Alhammoud  et  ah,  2005). 

The  LB  communicates  with  the  Ionian  Basin  through  the  Cretan  passage 
which  is  300  km  wide  and  more  than  2000  m  deep,  on  average,  between 
the  Cretan  and  Libyan  coasts.  To  the  northwest,  it  is  connected  to  the 
Aegean  Sea  through  three  passages  between  the  islands  of  Crete  and 
Karpathos  (Kassos  Strait,  width  67  km,  sill  depth  350  m),  the  islands  of 
Karpathos  and  Rhodes  (Karpathos  Strait,  43  km  and  550  m),  and  Rhodes 
and  Turkey  (17  km  and  350  m).  The  southeastern  continental  shelf  at  the 
Nile  delta  is  separated  into  two  shallow  areas:  Mersa-Matruh  to  the  west 
and  Shikmona  to  the  east  south  of  Cyprus.  (Alhammoud  et  ah,  2005) 


Eigure  4.  Bottom  Topography  of  the  Eastern  Mediterranean  Sea 


15  20  25  30  35 


Source:  Alhammoud,  B.,  K.  Beranger,  L.  Mortier,  M.  Crepon,  and  I.  Dekeyser,  2005: 
Surface  circulation  of  the  Levantine  Basin:  comparison  of  model  results  with 
observations.  Prog.Oceanogr.,  66,  299-320. 
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B,  WATER  MASSES 

The  processes  of  comprehensive  pertinence  for  ocean  climate  motions  cover 
features  of  the  thermohaline  circulations,  distribution,  water  mass  formations,  mixing  and 
transformation.  The  Atlantic  origin  fresh  water  comes  from  the  Atlantic  Ocean  to  the 
Mediterranean  Sea  through  the  Strait  of  Sicily.  Due  to  the  wind  and  other  forces,  this 
water  mass  spreads  and  evaporation  surpasses  precipitation  at  the  surface.  So,  the  water 
mass  evolves  into  much  saline  and  dense  water  mass.  The  cooling  develops  in  winter 
season  and  heating  develops  in  summer  season.  During  storm  developments,  especially  in 
March,  Levantine  Intermediate  Water  (LIW)  forms  in  the  Levantine  basin.  LIW  flows 
and  circulates  from  surface  to  a  few  hundred  meters  depth.  Some  of  LIW  masses  disperse 
inside  the  basin  and  some  leave  underneath  the  Atlantic  Water  (AW)  through  the  Strait  of 
Sicily  (Robinson  et  ah,  1992). 

LIW  is  the  most  important  origin  of  the  Mediterranean  water  masses  in  the 
Atlantic  (Robinson  et  ah,  1992).  The  representation  of  flow  through  the  Strait  of  Gibraltar 
is  seen  in  Figure  5. 


Figure  5.  Flow  Representation  through  the  Strait  of  Gibraltar 


Depth 

(m) 


Evaporation 


Atlantic  Ocean 


Mediterranean 


Source:  The  National  Oceanography  Centre  Southampton,  2015:  Mediterranean  Outflow 
Sill.  Accessed  04  October  2015.  [Available  online  at  http://www.euroargo-edu.org/img/ 
medoutflsill.png] 
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The  Mediterranean  Sea’s  thermohaline  cireulation  is  mostly  deseribed  by  an  open 
vertical  cell  covering  the  whole  basin  and  two  closed  vertical  cells  confined  by  the 
eastern  and  western  sub-basins  (Lascaratos  et  ah,  1999).  The  open  thermohaline  cell 
represents  the  shift  of  the  AW,  which  is  also  known  as  modified  Atlantic  Water  (MAW), 
to  the  LIW.  MAW  spreads  in  the  surface  layer,  while  the  LIW  circulates  in  the 
intermediate  layer.  The  two  closed  vertical  cells  represent  the  shift  of  intermediate  and 
surface  water  to  the  deep  waters,  which  are  Eastern  Mediterranean  Deep  Water  (EMDW) 
and  Western  Mediterranean  Deep  Water  (WMDW)  (Lascaratos  et  ah,  1999).  The 
thermohaline  circulation  sketch  of  the  Mediterranean  Sea  is  seen  in  Eigure  6. 


Eigure  6.  Thermohaline  Circulation  of  the  Mediterranean  Sea 


Formation  of  Western 
Mediterranean  Deep  Water 
(WMDW) 


Formation  of  Eastern 
Mediterranean  Deep 


Transport  of  Levantine 
Intermediate  Waters  (LIW) 


Source:  Demirov,  E.  K.,  and  N.  Pinardi,  2007:  On  the  relationship  between  the  water 
mass  pathways  and  eddy  variability  in  the  Western  Mediterranean  Sea.  Journal  of 
Geophysical  Research,  112,  1-21,  doi:10.1029/2005JC003174. 
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Three  main  permanent  water  masses  are  observed  in  the  Eastern  Mediterranean 
Sea.  They  are:  modified  Atlantic  Water  (MAW),  Levantine  Intermediate  Water  (LIW), 
and  Eastern  Mediterranean  Deep  Water  (EMDW)  (Alhammoud  et  ah,  2005). 

1.  Modified  Atlantic  Water  (MAW) 

Modified  Atlantic  Water  (MAW)  spreads  towards  the  Eastern  Mediterranean 
from  the  Strait  of  Sicily  in  the  surface  layer.  It  adjusts  the  mass  loss  of  the  Eastern 
Mediterranean  Deep  Water  by  cause  of  evaporation.  As  it  spreads  towards  the  Eastern 
Mediterranean,  its  salinity  is  intensified  and  its  depth  range  is  declined  due  to  the  high 
evaporation  (Alhammoud  et  ah,  2005). 

The  Modified  Atlantic  Water  is  comparably  less  saline  and  warmer.  The  MAW 
primarily  affects  the  upper  thermocline,  which  is  0-1 50m  (Manca  et  ah,  2004).  The 
MAW  is  characterized  as  an  extensive  homogeneous  layer  and  its  signature  is  seasonal. 
The  minimum  salinity  of  the  MAW  is  nearly  37.2  psu,  which  is  spotted  at  around  50 
meters  during  summer  and  its  salinity  minimum  is  located  near  the  surface  during  the 
winter.  The  temperature  range  is  between  15  °C  and  17  °C  and  the  salinity  range  is 
between  37.2  psu  and  37.8  psu  up  to  100  meters  depth  (Drago  et  ah,  2010). 

2,  Levantine  Intermediate  Water  (LIW) 

Levantine  Intermediate  Water  (LIW)  is  the  most  crucial  water  mass  in  the  Eastern 
Mediterranean  (Alhammoud  et  ah,  2005).  It  influences  not  only  the  whole  Mediterranean 
Sea  but  also  the  Atlantic  Ocean.  It  is  a  homogeneous  water  mass  below  the  Atlantic 
Water  and  circulates  nearly  between  300  and  700  meters.  LIW  occurs  in  the  Levantine 
basin  in  winter  season  due  to  the  transformation  of  the  Atlantic  Water.  As  LIW  spreads 
westward,  its  salinity  declines  because  of  the  mixing  with  the  bordering  water  masses 
(Alhammoud  et  ah,  2005). 

LIW  is  comparably  much  saline  and  colder.  Its  salinity  values  are  between  38.73 
psu  and  38.78  psu,  and  temperature  values  are  between  13.75  °C  and  13.92  °C  at  a  mean 
depth  of  nearly  300  meters  at  the  Strait  of  Sicily.  Its  salinity  maximum  occurs  in  the 
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western  and  south-western  of  Malta.  The  recurrenee  time  of  LIW  in  the  Strait  of  Sicily  is 
predicted  as  9  months  (Drago  et  ah,  2010). 

3,  Eastern  Mediterranean  Deep  Water  (EMDW) 

Eastern  Mediterranean  Deep  Water  circulates  nearly  below  800  meters  depth.  It  is 
dense  water  mass  and  relatively  less  saline  and  colder.  It  forms  in  the  northern  section  of 
the  eastern  basin  due  to  the  strong  cooling.  Its  mean  salinity  value  is  around  38.7  psu  and 
its  mean  temperature  value  is  nearly  13.6  °C  (Alhammoud  et  ah,  2005).  Vertical 
distribution  of  the  water  masses  in  the  Mediterranean  Sea  is  seen  in  Figure  7. 


Figure  7.  Vertical  Distribution  of  Water  Masses  in  the  Mediterranean 

Sea 


Depth 


Mediterranean  Sea  water  masses:  vertical  distribution 

Strdits  of  Gibraltar  Straits  of  Skily 


u 

Modified 

, 

100 

36.2 

Atlantic  Water 
(MAW) 

37.0 

Levantine 

200 

,  13.0 

38.4 

Levantirse 

Intermediate 

Water 

,  14.2 

3&8 

Intermediate 

Water 

(LIW) 

300 

(LIW) 

38.6 


>5.5 

39.1 


400 

SOO 

1000 

2000 

3000 

4000 


West 

12/  Mediterranean 
.(^.1  DeepWatei 
:wMOWi 


Water  movement 


14J  IempetQture{*C) 
S8.8  SatinitvirsU) 


lOMV 


I 

0* 


10*£ 


East 

Mediterrarsean 
Deep  Water 
(EMOW' 


I 

20^ 


;3.6 

38? 


Notr  DrpdilinbnottOujte 

r%  Ptm  tK  jl  Unk. 

Sewc*  adipM horn M.andMetor.G 
U  A  Munwncil  Slwt)r  ol  itw  MedMnanMn  Sm 
OKUMioa  AMncan  HeMofologKJt  Soaety.  I  MS 

T 

30*E  Longitude 


Source:  GRID-Arendal,  2013:  Mediterranean  Sea  water  masses:  vertical  distribution. 
Accessed  08  October  2015.  [Available  online  at  http://www.grida.no/graphicslib/detail/ 
mediterranean-sea-water-masses- vertical-distribution_d84b] 
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C.  WIND  FIELDS 

The  Mediterranean  basin  is  nearly  bounded  by  mountains.  The  meteorologieal 
events  in  this  area  are  affeeted  by  mountain  ehains  and  defined  by  the  other  sub-basins. 
The  importance  of  the  sub-basins  comes  from  the  synoptic -sc ale  airflow,  which  usually 
communicates  with  the  surrounding  orography  (Zecchetto  and  De  Biasio,  2007).  Spatial 
distribution  of  the  annual  mean  wind  speed  and  wind  direction  of  the  Mediterranean  Sea 
is  seen  in  Figure  8. 


Figure  8.  Spatial  Distribution  of  the  Annual  Mean  Wind  Speed  (Upper)  and  Wind 
Direction  (Lower)  of  the  Mediterranean  Sea 

Spatial  distribution  of  mean  value  of  .  Annual 


Spatial  distribution  of  wind  directionality.  Annual 


Source:  Gaillard,  P.,  and  Coauthors,  2004:  Wind  and  wave  atlas  of  the  Mediterranean 
Sea.  Accessed  15  October  2015.  [Available  online  at  http://users.ntua.gr/mathan/pdl7 
Pages_from%20_WIND_WAVE_ATLAS_MEDITERRANEAN_SEA_2004.pdf|. 
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The  regional  prevailing  winds  in  the  Mediterranean  region  are  westerlies  with  a 
northwesterly  tilt  and  cross  over  the  Mediterranean  Sea  longitudinally.  The 
Mediterranean  region  is  characterized  by  the  mountain  ranges  and  this  also  assists  to 
describe  the  wind  climate  of  the  basin  and  its  sub-basins.  During  the  winter  months, 
cyclonic  storms  are  leading  systems  in  the  Mediterranean  region  (Heidom,  2007). 

Mediterranean  region  is  the  center  place  for  the  cyclogenesis,  beginning  of  low 
pressure  systems  during  the  winter  season.  It  is  believed  that  the  highest  consolidation  of 
regionally  grown  cyclones  on  the  globe  occurs  in  the  Mediterranean  region,  partly  during 
the  winter  season  (Heidom,  2007).  Infrequently,  some  of  these  cyclones  grow  into  giant 
storms  looking  like  Atlantic  hurricanes;  however,  they  are  not  actual  tropical  storms. 
Although  the  Mediterranean  cyclones  are  tinier  and  mostly  less  damaging  than  the 
Atlantic  and  Pacific  tropical  storms,  these  Mediterranean  cyclones  are  very  much  alike  to 
the  tropical  cyclones  and  polar  lows  in  terms  of  the  properties  in  the  stmcture.  Some  of 
the  important  Mediterranean  cyclones  developed  into  hurricane-like  conditions  in  the 
region  during  the  last  sixty  years  (Heidom,  2007).  The  Mediterranean  Storm,  which 
occurred  in  the  Mediterranean  region  in  15  January  1995,  is  seen  in  Figure  9. 

Figure  9.  Hurricane-like  Cyclone  System  Occurred  in  15  January 
1995  in  the  Mediterranean  Region 


Mediterranean  Storm 
15  January  1995 


Source:  Heidom,  K.  C.,  2007:  The  Mediterranean:  Birthplace  of  the  Winds.  Accessed  24 
October  2015.  [Available  online  at  http://www.islandnet.com/~see/weather/almanac/ 
arc2007/alm07mar.htm] 
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The  general  wind  direetionality  in  the  Mediterranean  Sea  is  from  north  and  north¬ 
west  as  seen  in  Figures  8,10,1 1,12  and  13,  respeetively. 

The  oeean  waves  oecur  on  the  region  from  north-west  to  south-east  due  to  winds 
in  the  Mediterranean  region.  Thus,  large  waves  are  observed  at  the  north  shores  of 
Algeria,  Tunisia,  and  Egypt.  These  areas  are  needed  to  be  surveyed  earefully  for  the 
safety  eoneems  of  the  coastal  structures  (Ancient  Ports,  2015). 

During  the  winter,  winds  are  generally  westerlies  and  north  westerlies  with 
average  velocities  of  12  to  17  knots.  While  the  south-west  winds  predominate  in  the 
Adriatic  Sea,  both  east  and  west  winds  are  common  between  the  Strait  of  Gibraltar  and 
Sardinia.  Generally,  the  north  and  north-east  winds  are  observed  around  the  Cyprus  and 
the  northeastern  Mediterranean  Sea.  The  Aegean  Sea  departs  from  the  general  pattern  in 
the  Eastern  Mediterranean  Sea,  and  north  winds  dominate  in  all  seasons  whereas  both 
north  and  sound  winds  are  recurrent  in  the  winter  (National  Geospatial-Intelligence 
Agency,  2011).  Spatial  distribution  of  the  winter  season  wind  direction  of  the 
Mediterranean  Sea  is  seen  in  Eigure  10. 


Eigure  10.  Spatial  Distribution  of  Winter  Season  Wind  Direction  of 

the  Mediterranean  Sea 

Spatial  distribution  of  wind  directionality.  Winter 


Source:  Gaillard,  P.,  and  Coauthors,  2004:  Wind  and  wave  atlas  of  the  Mediterranean 
Sea.  Accessed  15  October  2015.  [Available  online  at  http://users.ntua.gr/mathan/pdf7 
Pages_from%20_WIND_WAVE_ATLAS_MEDITERRANEAN_SEA_2004.pdf|. 
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Apart  from  the  northeastern  Mediterranean  Sea,  the  general  wind  pattern  has 
slight  variation  in  the  spring  season,  where  south-west  winds  are  noticeable  at  this  time 
with  average  wind  velocities  9  to  13  knots  in  the  spring  season.  Spatial  distribution  of  the 
spring  season  wind  direction  of  the  Mediterranean  Sea  is  seen  in  Figure  1 1 . 


Figure  1 1 .  Spatial  Distribution  of  Spring  Season  Wind  Direction  of  the 

Mediterranean  Sea 

Spatial  distribution  of  wind  directionality.  Spring 


Source:  Gaillard,  P.,  and  Coauthors,  2004:  Wind  and  wave  atlas  of  the  Mediterranean 
Sea.  Accessed  15  October  2015.  [Available  online  at  http://users.ntua.gr/mathan/pdl7 
Pages_from%20_WIND_WAVE_ATLAS_MEDITERRANEAN_SEA_2004.pdf|. 

The  average  wind  velocities  are  5  to  9  knots  over  the  area  in  the  summer  season. 
The  general  wind  pattern  stays  unchanging  other  than  the  Aegean  Sea,  where  south  winds 
are  infrequent  (National  Geospatial-Intelligence  Agency,  2011).  Fall  season  shows  little 
difference  in  the  general  wind  pattern  excluding  the  southeastern  Mediterranean  Sea, 
where  winds  are  fluctuating.  The  average  wind  velocities  grow  to  7  to  10  knots  in  this 
season  (National  Geospatial-Intelligence  Agency,  2011).  Spatial  distributions  of  the  wind 
directionality  for  summer  and  autumn  seasons  are  depicted  in  Figure  12  and  Figure  13, 
respectively. 

Arrows,  which  represent  the  direction  of  the  winds,  only  resemble  the  winds 

larger  than  6  meter  per  second  and  frequent  larger  than  10%.  The  arrow  lengths  are 

relative  to  the  occurrence  frequency  (Gaillard  et  ah,  2004). 
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Figure  12.  Spatial  Distribution  of  Summer  Season  Wind  Direction  of 

the  Mediterranean  Sea 


Spatial  distribution  of  wind  directionality.  Summer 
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Source:  Gaillard,  P.,  and  Coauthors,  2004:  Wind  and  wave  atlas  of  the  Mediterranean 
Sea.  Accessed  15  October  2015.  [Available  online  at  http://users.ntua.gr/mathan/pdf7 
Pages_from%20_WIND_WAVE_ATLAS_MEDITERRANEAN_SEA_2004.pdf|. 


Figure  13.  Spatial  Distribution  of  Autumn  Season  Wind  Direction  of 

the  Mediterranean  Sea 


Spatial  distribution  of  wind  directionality.  Autumn 


Source:  Gaillard,  P.,  and  Coauthors,  2004:  Wind  and  wave  atlas  of  the  Mediterranean 
Sea.  Accessed  15  October  2015.  [Available  online  at  http://users.ntua.gr/mathan/pdf7 
Pages_from%20_WIND_WAVE_ATLAS_MEDITERRANEAN_SEA_2004.pdf|. 
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The  differential  heating  of  the  land  and  sea,  whieh  is  observed  in  the  warm  season 
and  oecasionally  apparent  in  the  eool  season,  leads  to  day-to-day  shift  of  the  land  and  sea 
breezes.  The  typical  sea  breezes  dominate  from  April  to  October,  onset  at  0700  or  0800, 
extending  to  a  maximum  around  1300  or  1400,  and  lasting  until  around  1800  or  1900. 
The  sea  breeze  appears  later  in  the  morning  in  the  spring  and  autumn;  while  its  beginning 
may  be  postponed  until  noon  in  the  winter  season.  The  land  and  sea  breezes  extend 
around  10  to  20  miles  from  the  coast.  While  the  mean  magnitudes  of  the  land  breezes  are 
5  to  9  knots,  the  sea  breezes  are  stronger,  reaching  to  1 1  to  16  knots  (National  Geospatial- 
Intelligence  Agency,  2011). 

1.  Etesian  Winds 

The  most  widespread  regional  wind  in  the  Levantine  basin  of  the  Eastern 
Mediterranean  Sea  is  the  Etesian  winds.  The  name  comes  from  the  Greek  origin  “etos,” 
which  meaning  “yearly”  to  indicate  the  Etesian  winds’  annual  endurance  (Zecchetto  and 
De  Biasio,  2007). 

The  regular  directions  of  the  Etesian  winds  are  from  north  and  northwestern.  The 
cold  air  and  clear  skies  are  observed  as  a  result  effect  of  Etesian  winds  between  the  late 
May  and  early  October.  Although  they  are  generally  dry,  cool  and  moderate;  they  may 
turn  into  a  windstorm  over  off-shore  areas,  extending  maximum  wind  speed  early 
afternoon  (WeatherOnline,  2015). 

The  high  pressure  system  over  the  Balkans  and  southeastern  Europe,  low  pressure 
system  over  Turkey,  and  evolution  of  the  south  Asian  warmth  over  the  Anatolian  plateau 
causes  the  occurrence  of  the  Etesian  winds  (WeatherOnline,  2015).  They  are  a 
component  of  the  Asian  monsoon  system.  Thermal  lows  deepening  over  Turkey  may 
cause  windstorm  type  Etesian  winds.  They  can  generate  important  seas  along  the 
Egyptian  and  Israeli  shores  due  to  the  lengthened  span  of  these  winds  (WeatherOnline, 
2015).  Eocation  of  Etesian  winds  is  seen  in  Eigure  14. 


17 


Figure  14.  Location  of  the  Etesian  Winds 


Source:  WeatherOnline,  2015:  Etesian  Winds.  Accessed  12  November  2015.  [Available 
online  at  http://www.weatheronline.co.uk/reports/wind/The-Etesian-Winds.htm] 


2,  Khamsin  Winds 

The  Khamsin  winds  develop  in  north  of  Africa  and  around  the  Eastern 
Mediterranean  region  periodically  between  the  end  of  winter  and  the  beginning  of 
summer.  They  are  overwhelming,  warm,  dry,  and  sandy  south  /  south-east  winds.  The 
name  “Khamsin”  comes  from  the  Arabic  originated  “Khamsun”  or  “Hamsin,”  which 
means  “fifty”  referring  to  the  duration  time  of  the  wind.  The  Khamsin  winds  appear 
advance  of  depressions,  which  maneuver  to  eastward  or  north-eastward  in  the 
Mediterranean  Sea  or  over  north  of  Africa  (WeatherOnline,  2015).  Settlement  of  the 
Khamsin  winds  in  the  Eastern  Mediterranean  Sea  is  seen  in  Eigure  15. 
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Figure  15.  Settlement  of  the  Khamsin  Winds  in  the  Eastern 
Mediterranean  Sea 


Source:  Met  Office,  2015:  Winds  of  the  world.  Accessed  12  November  2015.  [Available 
online  at  http://www.metoffice.gov.uk/learning/wind/wind-names] 


D.  NUTRIENTS 

The  Mediterranean  Sea  is  generally  distinguished  by  low  nutrient  concentration, 
especially  in  the  Eastern  Mediterranean  (Tanhua  et  ah,  2013).  The  surface  layers  are 
mostly  quite  nutrient  exhausted.  Due  to  low  nutrient  concentration,  Mediterranean  can  be 
classified  as  oligotrophic  or  indeed  ultra-oligotrophic  sea.  The  nutrient  exhausted  surface 
layer  of  the  Eastern  Mediterranean  is  thicker  and  the  deeper  comparing  to  the  Western 
Mediterranean  (Tanhua  et  ah,  2013). 

Deep  water  with  greater  nutrient  values  flows  towards  the  Atlantic  Ocean,  while 
the  Atlantic  Water  with  lower  nutrient  values  flows  towards  the  Mediterranean  Sea.  This 
circulation  pattern,  known  as  anti-estuarine  circulation,  is  one  of  the  main  reasons  for 
lower  nutrient  concentration  in  the  Mediterranean  Sea  (Tanhua  et  ah,  2013). 

The  apparent  oxygen  utilization  (AOU)  rate  is  nearly  70-80  ^8  .  Even  at 

the  oxygen  minimum  layer  of  the  intermediate  depths,  the  oxygen  concentration  rate  is 

nearly  IgOE^^^^^S  .  Considering  this  rate  difference,  the  Mediterranean  is  known  as  a 


19 


well  oxygenated  sea.  This  high  oxygen  eoncentration  in  the  Mediterranean  Sea  is  due  to 
the  dynamic  and  quick  circulation  of  the  Mediterranean  Sea  (Tanhua  et  ah,  2013). 

The  dissolved  oxygen  and  nutrient  concentrations  in  the  Eastern  Mediterranean 
Sea  show  largely  homogenous  distribution  in  the  upper  layer.  The  surface  waters  are 
almost  saturated  with  the  dissolved  oxygen.  The  distribution  of  nutrients  in  the  Levantine 
basin  is  seen  in  Figure  16  (Krom  et  ah,  2014). 


Figure  16.  Vertical  Sections  of  Dissolved  Oxygen  and  Nutrients  along 
West-East  Cross  Section  of  the  Levantine  Basin 
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Source:  Krom,  M.,  N.  Kress,  I.  Berman-Frank,  and  E.  Rahav,  2014:  Past,  Present  and 
Future  Patterns  in  the  Nutrient  Chemistry  of  the  Eastern  Mediterranean.  The 
Mediterranean  Sea  Its  history  and  present  challenges,  S.  Goffredo,  and  Z.  Dubinsky, 
Eds.,  Springer,  49-68. 


The  distribution  of  the  dissolved  oxygen  concentrations  is  large  at  the  surface.  As 
depth  increases,  the  dissolved  oxygen  concentration  reaches  to  the  slightest  values  at  the 
500-1 500m  layer.  Then,  it  rises  towards  the  bottom  of  the  Levantine  basin,  or  down  to 
2500m  (Krom  et  ah,  2014).  Composite  depth  profdes  of  dissolved  oxygen  and  nutrients 
in  the  Levantine  basin  are  seen  in  Figure  17. 
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The  nitrate  depth  profile  is  nearly  opposite  of  the  dissolved  oxygen  profile.  The 
nitrate  eoneentration  is  small  at  the  surface.  It  reaches  to  the  greatest  values  at  the  500- 
1500m  layer,  and  then  it  reduces  towards  the  bottom  of  the  Levantine  basin,  or  down  to 
2500m.  The  silicic  acid  depth  profile  shows  a  similar  pattern  with  the  nitrate  depth  profile 
(Krom  et  ah,  2014). 


Figure  17.  Composite  Depth  Profiles  of  Dissolved  Oxygen  and 
Nutrients  in  the  Levantine  Basin 
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Source:  Krom,  M.,  N.  Kress,  I.  Berman-Frank,  and  E.  Rahav,  2014:  Past,  Present  and 
Future  Patterns  in  the  Nutrient  Chemistry  of  the  Eastern  Mediterranean.  The 
Mediterranean  Sea  Its  history  and  present  challenges,  S.  Goffredo,  and  Z.  Dubinsky, 
Eds.,  Springer,  49-68. 
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III.  DATA 


A,  HISTORICAL  BACKGROUND  OF  THE  DATASETS 

The  data  used  in  this  study  were  obtained  from  the  World  Oeean  Database  2013 
(WOD13),  whieh  was  developed  by  the  Oeean  Climate  Laboratory  (OCL)  at  the  National 
Oceanographie  Data  Center  (NODC),  Silver  Spring,  Maryland,  USA.  The  scientifically 
quality-controlled  in-situ  sub-surface  and  surface  measurements  were  included  in  the 
WOD13  (Johnson  et  ah,  2013). 

The  first  release  of  these  dataset  series  was  the  World  Ocean  Atlas  1994 
(WOA94).  It  was  comprised  of  temperature,  salinity,  oxygen,  phosphate,  nitrate,  and 
silicate  profiles.  The  second  release  of  the  dataset  was  World  Ocean  Database  1998 
(WOD98).  This  release  renewed  the  previous  one  by  adding  new  data  as  well  as  new 
variables.  The  third  release  of  the  dataset  was  World  Ocean  Database  2001(WOD01). 
The  data,  which  were  collected  from  new  instruments  types  such  as  Undulating  Ocean 
Recorders  and  profiling  floats,  were  covered  in  WODOl. 

The  fourth  release  of  the  dataset  was  World  Ocean  Database  2005  (WOD05). 
This  dataset  updated  the  previous  one  by  containing  data  from  new  variables  and  from 
one  new  instrument  that  is  gliders.  The  fifth  update  of  the  dataset  was  World  Ocean 
Database  2009  (WOD09).  It  added  a  new  feature  that  was  quality  control.  The  data  in  this 
release  were  updated  every  three  months  (Boyer  et  ah,  2013). 

The  newest  release  of  these  series  is  World  Ocean  Database  2013  (WOD13).  In 
this  release,  the  standard  depth  levels  were  expanded  from  40  to  138  to  have  better 
resolution.  The  comparison  of  amount  of  the  data  in  WOD13  with  the  previous  versions 
of  World  Ocean  Database  (WOD)  and  standard  depth  levels  are  seen  in  Table  2  and 
Table  3,  respectively  (Boyer  et  ah,  2013). 
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Table  2.  Comparison  of  amount  of  Data  in  WOD13  with  Previous  Versions 

ofWOD 


Dataset 

NODC 

(1974)’ 

NODC 

(1991)* 

WOA94 

WOD98 

WOD01 

WOD05 

WOD09 

WOD13 

OSD  ^ 

425,000 

783,912 

1,194.407 

1,373,440 

2,121,042 

2,258,437 

2,541,298 

3,115,552 

ctd' 

na 

66,450 

89,000 

189,555 

311,943 

443,953 

641,845 

848,911 

mbt' 

775,000 

980,377 

1,922,170 

2,077,200 

2,376,206 

2,421,940 

2,426,749 

2,425,607 

XBT 

290,000 

704,424 

1,281,942 

1,537,203 

1,743,590 

1,930,413 

2,104.490 

2,211,689 

MRB 

na 

na 

na 

107,715 

297,936 

445,371 

566.544 

1.411,762 

DRB 

na 

na 

na 

na 

50,549 

108,564 

121,828 

251,712 

PFL 

na 

na 

na 

na 

22,637 

168,988 

547,985 

1,020,216 

UOR 

na 

na 

na 

na 

37,645 

46.699 

88,190 

88,190 

APB 

na 

na 

na 

na 

75,665 

75,665 

88.583 

1,713,132 

GLD 

na 

na 

na 

na 

na 

338 

5,857 

103,798 

Total 

Stations 

1,490.000 

2,535.163 

4,487,519 

5.285.113 

7,037,213 

7.900,368 

9,155.099 

13,190,569 

Plankton 

na 

na 

na 

83,650 

142,900 

150,250 

218,695 

242,727 

sur‘ 

na 

na 

na 

4,743 

9,178 

9,178 

9,289 

Source:  Johnson,  D.,  T.  Boyer,  H.  Garcia,  R.  Locamini,  O.  Baranova,  and  M.  Zweng, 
2013:  World  Ocean  Database  2013  User’s  Manual  NODC  Internal  Report  22,  163  pp, 
[Available  online  at  http://www.nodc.noaa.gov/OC5AVOD13/docwodl3.html] 


Table  3.  Standard  Levels  and  Depths  in  Meters 


Depth 

Level # 

Depth 

Level  ** 

Depth 

Level # 

Depth 

Level # 

O 

1 

475 

36 

2400 

71 

5900 

106 

5 

2 

500 

37 

2500 

72 

6000 

107 

lO 

3 

550 

38 

2600 

73 

6100 

108 

15 

4 

600 

39 

2700 

74 

6200 

109 

20 

5 

650 

40 

2800 

75 

6300 

no 

25 

6 

700 

41 

2900 

76 

6400 

in 

30 

7 

750 

42 

3000 

77 

6500 

112 

35 

8 

800 

43 

3100 

78 

6600 

113 

40 

9 

850 

44 

3200 

79 

6700 

114 

45 

lO 

900 

45 

3300 

80 

6800 

115 

50 

11 

950 

46 

3400 

81 

6900 

116 

55 

12 

lOOO 

47 

3500 

82 

7000 

117 

60 

13 

1050 

48 

3600 

83 

7100 

118 

65 

14 

IlOO 

49 

3700 

84 

7200 

119 

70 

15 

1150 

50 

3800 

85 

7300 

120 

75 

16 

1200 

51 

3900 

86 

7400 

121 

80 

17 

1250 

52 

4000 

87 

7500 

122 

85 

18 

1300 

53 

4100 

88 

7600 

123 

90 

19 

1350 

54 

4200 

89 

7700 

124 

95 

20 

1400 

55 

4300 

90 

7800 

125 

lOO 

21 

1450 

56 

4400 

91 

7900 

126 

125 

22 

1500 

57 

4500 

92 

8000 

127 

150 

23 

1550 

58 

4600 

93 

8100 

128 

175 

24 

1600 

59 

4700 

94 

8200 

129 

200 

25 

1650 

60 

4800 

95 

8300 

130 

225 

26 

1700 

61 

4900 

96 

8400 

131 

250 

27 

1750 

62 

5000 

97 

8500 

132 

275 

28 

1800 

63 

5100 

98 

8600 

133 

300 

29 

1850 

64 

5200 

99 

8700 

134 

325 

30 

1900 

65 

5300 

lOO 

8800 

135 

350 

31 

1950 

66 

5400 

lOl 

8900 

136 

375 

32 

2000 

67 

5500 

102 

9000 

137 

400 

33 

2100 

68 

5600 

103 

425 

34 

2200 

69 

5700 

104 

450 

35 

2300 

70 

5800 

105 

Source:  Johnson,  D.,  T.  Boyer,  H.  Garcia,  R.  Locamini,  O.  Baranova,  and  M.  Zweng, 
2013:  World  Ocean  Database  2013  User’s  Manual  NODC  Internal  Report  22,  163  pp, 
[Available  online  at  http://www.nodc.noaa.gov/OC5/WOD13/docwodl3.html] 
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B, 


THE  FEATURES  OF  THE  DATASETS  IN  WOD13 


The  data  in  WOD13  are  comprised  of  eleven  datasets,  each  of  which  is 
characterized  by  the  instruments  that  collect  the  data.  Moreover,  the  variables  in  the 
WOD13  database  were  measured  both  at  observed  depth  levels  and  at  138  standard  depth 
levels.  The  datasets  in  WOD13  are  listed  in  Table  4. 


Table  4.  Datasets  in  the  WOD13 


DATASETS 

DATASETS  INCLUDES 

OSD 

Ocean  Station  Data,  Low-resolution  CTD/XCTD,  Plankton  data 

CTD 

High-resolution  Conductivity-Temperature-Depth  /  XCTD  data 

MBT 

Mechanical  /  Digital  /  Micro  Bathythemiograph  data 

XBT 

Expendable  Bathythermograph  data 

SUR 

Surface-only  data 

APB 

Autonomous  Pinniped  data 

MRB 

Moored  buoy  data 

PEL 

Profiling  float  data 

DRB 

Drifting  buoy  data 

UOR 

Undulating  Oceanographic  Recorder  data 

GLD 

Glider  data 

Source:  Johnson,  D.,  T.  Boyer,  H.  Garcia,  R.  Locarnini,  O.  Baranova,  and  M.  Zweng, 
2013:  World  Ocean  Database  2013  User’s  Manual  NODC  Internal  Report  22,  163  pp, 
[Available  online  at  http://www.nodc.noaa.gov/OC5AVOD13/docwodl3.html] 


1,  Ocean  Station  Dataset  (OSD) 

Ocean  Station  Data  includes  measurements  that  are  conducted  via  research 
vessels.  The  reversing  thermometers  are  used  to  measure  the  temperature  and  special 
bottles  are  used  to  gather  the  seawater  samples  (Boyer  et  ah,  2013).  The  salinity,  oxygen, 
nutrient  and  chlorophyll  are  the  variables  that  are  measured  via  collected  water  samples 
(Boyer  et  ah,  2013). 
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2,  Conductivity  Temperature  Depth  (CTD)  Dataset 

a.  High-Resolution  CTD  Data 

CTD  dataset  are  comprised  of  oxygen,  chlorophyll,  temperature  and  salinity 
profiles  (Boyer  et  ah,  2013).  If  the  depth  increment  of  the  data  is  less  than  two  meters,  it 
is  classified  as  high-resolution  CTD  (Boyer  et  ah,  2013). 

b.  Low-Resolution  CTD  Data 

The  pressure  sensor,  which  is  used  to  convert  the  measured  pressure  to  depth;  a 
temperature  device;  and  a  conductivity  sensor,  which  is  used  to  evaluate  the  salinity 
measurements,  is  the  combination  to  refer  the  Conductivity-Temperature -Depth  (CTD) 
instruments.  If  the  depth  increment  is  more  than  two  meters,  it  is  classified  as  low- 
resolution  CTD  (Johnson  et  ah,  2013). 

3,  Mechanical/Digital/Micro  Bathythermograph  (MBT)  Dataset 

The  MBT  data  provides  temperature  information  as  a  function  of  depth.  The  MBT 
data  is  collected  by  using  Mechanical  Bathythermograph  (MBT),  Digital 
Bathythermograph  (DBT)  and  Micro  Bathythermograph  (micro  BT)  instruments 
(Johnson  et  ah,  2013). 

4,  Expendable  Bathythermograph  (XBT)  Dataset 

The  deployment  of  the  XBTs  dates  back  to  the  1966.  The  expendable 
bathythermograph  is  generally  used  instead  of  MBT  in  most  measurements.  It  has  a 
thermistor,  which  measures  temperature  versus  depth.  This  instrument  can  be  submerged 
from  ship,  airplane,  or  submarine.  There  is  a  command  unit  for  recoding  the  temperature 
information  over  the  copper  wire  (Boyer  et  ah,  2013). 
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5. 


Profiling  Float  (PFL)  Dataset 


Temperature,  salinity,  and  pressure  data,  which  is  gathered  from  the  drifting 
profiling  floats  such  as  Profiling  Autonomous  Lagrangian  Circulation  Explorer,  free- 
drifting  hydrographic  profiler.  Sounding  Oceanographic  Lagrangian  Observer  and 
Autonomous  Profiling  Explorer,  is  included  in  PEL  dataset  (Johnson  et  al.,  2013).  The 
Argo  Project  is  the  main  contributor  to  the  PPL  dataset  in  the  WOD13  (Johnson  et  ah, 
2013). 


6,  Glider  (GLD)  Dataset 

Data,  from  reusable  autonomous  underwater  vehicles  (AUVs),  is  included  in  the 
Glider  dataset.  AUVs  are  designed  to  reach  to  certain  depths  and  measure  the 
temperature,  salinity,  current  and  other  quantities  (Boyer  et  ah,  2013).  Different  types  of 
gliders  are  seen  in  Pigure  18. 


Pigure  18.  Different  Types  of  Gliders 


Source:  Integrated  Ocean  Observing  System,  2014:  Gliders/Autonomous  Underwater 
Vehicles.  Accessed  15  November  2015.  [Available  online  at  http://www.ioos.noaa.gov/ 
glider/welcome.html] 
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C.  CHARACTERISTICS  OF  THE  DATA  USED  IN  THE  THESIS 

The  data  used  in  the  thesis  are  from  1960  to  2013.  The  data  obtained  from 
WOD13  are  re-organized  according  to  the  geographical  area  of  interest,  which  is  Eastern 
Mediterranean  Sea.  According  to  the  annual  distribution  of  temperature  profiles,  temporal 
distribution  of  temperature  in  the  Eastern  Mediterranean  Sea,  seen  in  Eigure  19,  shows 
that  the  maximum  number  of  temperature  profiles,  which  was  gathered  in  1968,  is  almost 
8,000;  while  the  minimum  number  of  the  temperature  profiles,  which  was  gathered  in 
1960,  is  almost  900.  According  to  the  monthly  distribution  of  temperature  profiles.  May 
and  October  have  the  maximum  number  of  temperature  profiles,  almost  18,000;  while 
December  has  the  minimum  number  of  the  temperature  profiles,  nearly  10,000. 


Eigure  19.  Temporal  Distribution  of  Temperature  Profiles  in  the 
Eastern  Mediterranean  Sea 
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It  is  clear  that  the  amount  of  temperature  profile  data  collected  along  the  Afriean 
eoast  is  less  compared  to  other  regions  in  the  Eastern  Mediterranean.  The  reliability  of 
the  study  depends  on  how  frequently  the  data  were  colleeted,  which  is  one  of  the  crucial 
aspects  in  determining  the  circulation  patterns.  The  analyses  ineluded  a  total  of  164,906 
temperature  profiles.  The  spatial  distribution  of  the  temperature  profiles  in  the  Eastern 
Mediterranean  region  is  seen  in  Eigure  20. 

Eigure  20.  Spatial  Distribution  of  Temperature  Profiles  in  the  Eastern 

Mediterranean 


According  to  the  temporal  distribution  of  the  salinity,  2005  has  the  maximum 
number  of  salinity  profiles,  nearly  3,700;  whereas  the  minimum  number  of  salinity 
profiles,  nearly  250,  was  obtained  in  1964.  The  temporal  distribution  histogram  of  the 
salinity  profiles  and  the  spatial  distribution  histogram  of  the  salinity  profiles  are  shown  in 
Figure  21  and  Figure  22,  respeetively.  October  has  the  maximum  number  of  salinity 
profiles,  nearly  6,000;  while  the  minimum  number  of  the  salinity  profiles,  nearly  2,900, 
was  produced  in  January.  The  analyses  ineluded  a  total  of  53,606  salinity  profiles. 
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Latitude 


Figure  21 .  Temporal  Distribution  of  Salinity  Profiles  in  the  Eastern 

Mediterranean 
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Figure  22.  Spatial  Distribution  of  Salinity  Profiles  in  the  Eastern 

Mediterranean 
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The  data  in  the  W0D13  were  eollected  via  many  oeeanographic  instruments  and 
WOD13  ineluded  many  datasets.  Types  of  the  datasets,  used  for  salinity  analyses,  are; 
Conductivity-Temperature -Depth  (CTD)  data.  Glider  (GLD)  data,  and  Ocean  Station 
Data  (OSD)  and  Profiling  Float  (PFL)  data.  The  number  of  the  monthly  and  annually 
salinity  profiles  and  their  distribution  according  to  the  different  type  of  datasets  are 
presented  in  Figure  23,  and  Figure  24,  respectively. 


Figure  23.  Monthly  Distribution  of  the  Different  Datasets  for  Salinity 
Profiles  in  the  Eastern  Mediterranean 
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Figure  24.  Annually  Distribution  of  the  Different  Datasets  for  Salinity 
Profiles  in  the  Eastern  Mediterranean 
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Types  of  the  datasets,  used  for  temperature  analyses,  are;  Conduetivity- 
Temperature-Depth  (CTD)  data.  Glider  (GLD)  data,  Oeean  Station  Data  (OSD),  Profiling 
Float  (PEL)  data.  Expendable  Bathythermograph  (XBT)  data  and  Meehanieal 
Bathythermograph  (MBT)  data.  The  number  of  the  monthly  and  annually  temperature 
profiles  and  their  distribution  according  to  the  different  type  of  datasets  are  presented  in 
Figure  25,  and  Figure  26,  respectively. 


32 


Frequency  Frequency  Frequency 


Figure  25.  Monthly  Distribution  of  the  Different  Datasets  for 
Temperature  Profiles  in  the  Eastern  Mediterranean 
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Figure  26.  Annually  Distribution  of  the  Different  Datasets  for 
Temperature  Profiles  in  the  Eastern  Mediterranean 
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IV.  METHODOLOGY 


A,  OPTIMAL  SPECTRAL  DECOMPOSITION  (OSD)  METHOD 

Data  assimilation  is  a  necessary  method  for  ocean  studies.  Moreover,  it  is  very 
significant  to  provide  the  realization  of  ocean  modeling  and  prediction.  The  Optimal 
Spectral  Decomposition  (OSD)  method  is  an  effective  approach  for  data  analysis.  The 
OSD  method  has  been  developed  for  ocean  data  analysis  to  establish  temporally  varying 
3-D  gridded  (T,  S)  fields  from  the  WOD13  profiles  (Chu  et  ah,  2015).  First  nine  basis 
functions  of  the  Eastern  Mediterranean  Sea  are  seen  in  Figure  27. 

The  basis  functions,  depending  only  on  the  topography  of  the  ocean  basin, 
are  the  eigenvectors  of  the  Laplacian  operator  with  the  same  lateral 
boundary  conditions  as  the  assimilated  variable  anomalies.  (Chu  et  ah, 

2015) 


Figure  27.  First  Nine  Basis  Functions  of  the  Eastern  Mediterranean 

Sea 
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Selection  of  basis  functions  |0A:(x,z)|  needs  to  satisfy  three  necessary 

conditions:  (i)  satisfaction  of  the  same  homogeneous  boundary  condition 
of  the  assimilated  variable  anomaly,  (ii)  orthonormal,  and  (iii) 
independence  of  the  assimilated  variable.  (Chu  et  al.,  2015) 

K  shows  the  number  of  basis  functions  and  N  shows  the  number  of  grid  points,  so 
orthonormal  can  be  expressed  as; 


N 


Y<3,{x,)0,{x,)  =  S„. 


(0.1) 


'kk' 


0  if  k^k' 
1  if  k  =  k' 


(0.2) 


Horizontal  grid  points  are  represented  by  xn  =  (xi,  yj),  n  =  1,  2,...,  N.  The 
following  matrix  represents  the  first  K  discrete  basis  functions  for  all  grid  points. 


0,(x,) 

0,(x,)  ■■ 

■  Z>,{x„) 

0Ax,) 

0;;  (.V;  ) 

■■ 

■■ 

(0.3) 
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B,  EMPIRICAL  ORTHOGONAL  FUNCTION  (EOF)  METHOD 


Empirical  Orthogonal  Function  (EOF)  method  has  been  used  in  meteorology  and 
oeeanography  in  order  to  interpret  the  spatial  and  temporal  analysis  related  to  the 
instability  of  variables  sueh  as  eurrents,  winds,  temperature  and  to  be  able  to  figure  out 
the  meehanisms  behind  those  variables  (Kaihatu  et  ah,  1998). 


The  empirical  orthogonal  eigenfunction  (EOF)  method,  first  introduced  to 
geophysics  by  Eorenz  (1956),  has  subsequently  been  used  to  enable 
analysis  of  data  with  eomplex  spatial/temporal  struetures.  Using  this 
method,  the  most  effieient  deeomposition  of  the  data  into  representative 
modes  is  determined  by  empirieally  finding  the  eigenfunctions  that  best 
deseribe  the  information.  (Kaihatu  et  ah,  1998) 

Retention  of  only  the  first  few  modes  may  eontain  a  signifieant  portion  of 
the  total  varianee,  leading  to  potentially  signifieant  storage  savings  if  not 
all  the  varianee  is  required.  Additionally,  each  mode  contains  phenomena 
with  differing  spatial/temporal  scales  and  thus  ean  be  isolated.  (Kaihatu  et 
ah, 1998) 

In  this  thesis,  EOF  method  allows  us  to  analyze  the  long  term  anomaly  of  the  heat 
eontent  and  freshwater  eontent.  Temperature  anomaly  field  T  is  obtained  by  subtraeting 
the  monthly  mean  T  from  temperature  data  T. 


f(x.,yj,z„Ti,t^)  =  T(x.,y., 


^k^^nL)-T[x„yj 


(0.4) 


f  is  modified  into  a  NxP  matrix  as  .  N  represents  the  number  of  grid 

points  and  P  represents  the  number  of  the  time  points.  In  this  study,  a  total  of  648  (54 
years  by  12  months)  time  points  were  used  for  P  value.  N  is  dependent  on  the  area  that  is 

used  for  analysis.  After  some  modifieations  of  the  equation  (1.4),  anomaly  field 

ean  be  written  as. 


=  (a U, ,  a )  (0.5) 

a 

where  temporal  variation  is  represented  byE*C^  (r,,^) ,  and  spatial  variation  is 
represented  by (aU;’  A )  • 
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Table  5  and  Table  6  show  the  varianees  of  first  four  EOF  modes  of  freshwater 
eontent  and  heat  eontent,  respeetively.  Aeeording  to  tables,  first  four  modes  represent  the 
nearly  50%  of  the  total  varianees,  whereas  the  remaining  26  modes  represent  the  other 
50%  of  the  total  varianee.  Thus,  it  is  elear  that  the  first  few  modes  are  much  noteworthy 
to  understand  the  general  distribution  of  variation. 

Table  5.  Freshwater  Content  Variances  of  First  Four  EOFs  Modes 


EOFs 

Variance  (%) 

Total  Variance  (%) 

1 

25.2 

25.2 

2 

11.9 

37.1 

3 

8.47 

45.57 

4 

6.3 

51.87 

Table  6.  Heat  Content  Variances  of  First  Four  EOFs  Modes 


EOFs 

Variance  (%) 

Total  Variance  (%) 

1 

27.8 

27.8 

2 

11 

38.8 

3 

6.02 

44.82 

4 

4.65 

49.47 

38 


C.  P-VECTOR  INVERSE  METHOD 

P-vector  inverse  method  is  used  to  obtain  the  absolute  geostrophic  veloeity  from 
temperature  and  salinity  (T,  S)  data.  There  are  three  required  physieal  basis  for  the  P- 
vector  inverse  method.  They  are:  mass  eonservation,  geostrophic  balance,  and  no  major 
cross-isopycnal  mixing  except  water  masses  in  contact  with  the  atmosphere  (Chu  et  ah, 
1998). 

The  conservation  of  mass  is  represented  by  the  equation, 

V-Vp=0  (0.6) 

where  p  represents  the  potential  density,  by  differentiating  this  equation  with  respect  to  z 
and  multiplying  with  the  coriolis  parameter  f,  conservation  of  potential  vorticity  equation 
can  be  produced  as, 

q  =  f^,V-Wq=0  (0.7) 

oz 

Equations  (1.6)  and  (1.7)  indicate  that  V  is  perpendicular  to  both  V/?  andV^ 
therefore,  the  velocity  V  is  parallel  to xVp .  (Chu  et  ah,  1998) 


So,  the  unit  vector  P  is  parallel  to  the  velocity  V ,  and  unit  vector  P  can  be 
produced  as. 


VpxVg 

VpxVq 


(0.8) 


As  a  result,  velocity  V  can  be  obtained  as. 


V  =  r[x,y,z)P  (0.9) 

where  r  represents  the  proportionality  factor  and  can  be  determined  by  set  of  algebraic 
equations  by  using  the  thermal  wind  relation  (Chu  et  ah,  1998). 
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V.  SEASONAL  VARIABILITY 


A,  TEMPERATURE  FIELD 

Figure  28  shows  the  total  mean  temperature  of  the  Eastern  Mediterranean  Sea  at 
different  depths  between  1960  and  2013.  The  sea  surfaee  temperature  changes  from  23  °C 
in  the  east  coast  of  the  Levantine  Sea  to  18  “C  in  the  entrance  of  the  Adriatic  Sea.  The 
temperature  values  have  slight  variation  from  surface  to  20m  depth.  Until  20m  depth,  the 
temperature  variation  is  nearly  1  -  2  °C.  The  mean  temperature  in  the  50m  depth  shows  a 
little  much  variation  comparing  to  20m  depth.  The  50m  depth  shows  transition  zone 
characteristics  between  the  surface  depth  and  the  150m  depth.  In  100m  depth,  the 
warmest  section  of  the  Eastern  Mediterranean  is  at  the  very  east  coast  of  the  Levantine 
Sea  with  a  temperature  of  almost  17  °C,  while  the  coldest  section  of  the  Eastern 
Mediterranean  is  at  the  Gulf  of  Taranto  and  Strait  of  Otranto  with  a  temperature  of  almost 
14.5  “C.  The  16  °C  isothermal  line  runs  from  Cretan  passage  to  Sfax,  passing  through  the 
north  coast  of  the  Egypt  and  Libya.  The  16  °C  and  15  °C  isothermal  lines  nearly  cover  all 
the  Central  Ionian  Sea.  The  1 50m  depth  is  the  beginning  of  the  intermediate  layer  of  the 
Eastern  Mediterranean  Sea.  Erom  150m  to  400m  depth,  the  temperature  difference 
between  the  east  part  of  the  Eastern  Mediterranean  and  the  west  part  of  the  Eastern 
Mediterranean  decreases  as  the  depth  goes  down.  In  400m  depth  layer,  temperature 
distribution  is  almost  uniform  through  the  basin  and  the  temperature  is  almost  15  °C. 
Seemingly,  there  is  no  another  isothermal  line  in  the  400m  depth  layer  except  the  very 
west  border  of  the  basin.  The  deep  layer  of  the  Eastern  Mediterranean  Sea  begins  at  the 
600m  depth.  There  is  nearly  no  temperature  change  at  600m  depth  layer  with  a 
temperature  of  almost  14  °C.  Below  600m,  temperature  remains  nearly  uniform  down  to 
the  bottom. 
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Figure  28.  Total  Mean  Temperature  (°C)  at  Different  Depths 
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Figure  29  shows  the  climatological  monthly  mean  temperature  at  the  surface. 
Temperature  drastically  drops  almost  3  °C  from  December  to  January.  As  time  goes  by, 
temperature  continues  to  decline  in  the  winter  season.  The  coldest  part  in  the  Eastern 
Mediterranean  is  the  entrance  of  the  Adriatic  with  a  temperature  of  13  °C.  Also,  February 
is  the  coldest  month  through  the  year.  In  the  spring  season,  March  draws  attention  with  a 
highly  similar  temperature  characteristic  of  February.  May  is  like  a  transition  month 
between  spring  and  summer  seasons.  Temperature  increases  to  19  °C  in  the  western  part 
and  to  22  °C  in  the  eastern  part  of  the  basin  in  May. 
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Relatively  eolder  part,  with  a  temperature  of  almost  19.5  °C,  is  the  entranee  of  the 
Strait  of  Gibraltar  for  June  in  summer  season.  Relatively  warmer  part,  with  a  temperature 
of  almost  24.5  “C,  is  the  eoast  of  Syria.  Temperature  increases  nearly  from  24  °C  to  27  °C 
throughout  the  whole  basin  in  July  and  increases  to  almost  29  °C  in  August.  August  is  the 
warmest  month  through  the  year.  East  part  of  the  basin  and  north  coasts  of  the  Africa  are 
the  warmest  areas  that  are  spotted  in  August.  With  the  beginning  of  the  fall  season, 
temperature  reduces  slowly  until  December. 

Figure  29.  Climatological  Monthly  Mean  Temperature  at  the  Surface 
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Figure  30  shows  the  elimatological  monthly  mean  temperature  at  the  150m  depth. 
The  mean  temperature  varies  almost  from  14.5  “C  to  17  °C  through  the  basin.  As 
observed  at  the  surface,  the  warmer  spots  are  detected  in  the  eastern  part  and  the  colder 
spots  are  detected  in  the  western  and  northern  part  of  the  basin. 


Figure  30.  Climatological  Monthly  Mean  Temperature  at  150m  Depth 
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Figure  3 1  shows  the  climatological  monthly  mean  temperature  at  the  600m  depth. 
The  temperature  distribution  is  almost  same  for  all  months.  Temperature  varies  from  13.5 
°C  to  14.5  °C.  While  the  eastern  part  of  the  basin  and  north  coasts  of  Africa  are  the 
warmer  spots  at  the  surface  and  150m  depth,  they  are  colder  spots  at  600m  depth.  Even 
though  a  significant  seasonal  variation  is  observed  at  the  surface  depth,  there  is  almost  no 
seasonal  temperature  variation  at  600m  depth. 

Figure  3 1 .  Climatological  Monthly  Mean  Temperature  at  600m  Depth 
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B,  SALINITY  FIELD 

Figure  32  shows  the  total  mean  salinity,  whieh  is  almost  uniformly  distributed 
through  the  basin  until  100m  depth  with  a  salinity  value  of  nearly  37  psu  in  the  entranee 
of  Strait  of  Gibraltar  and  with  a  salinity  value  of  nearly  39  psu  in  the  Levantine  Basin 
north-east  eomer,  where  Turkey  south  eoast  and  Syria  west  coast  intersects. 

Figure  32.  Total  Mean  Salinity  (psu)  at  Different  Depths 
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Between  100m  and  200m,  western  part  of  the  basin  around  Strait  of  Sicily  and 
Malta  becomes  much  saline  with  a  salinity  value  of  38  to  38.8  psu,  while  the  eastern  part 
of  the  basin  stays  same  with  a  salinity  value  of  nearly  39  psu.  From  200m  to  600m  depth, 
the  salinity  variation  is  very  slight  through  the  basin.  Both  eastern  and  western  parts  of 
the  basin  become  almost  equally  salted  with  a  mean  salinity  value  of  nearly  38.8  psu. 
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Below  600m  depth,  there  is  almost  no  variation  in  the  salinity  distribution  through  the 
basin. 

Figure  33  shows  the  climatological  monthly  mean  salinity  at  the  surface.  There  is 
no  significant  seasonal  salinity  variation  for  winter  and  spring  seasons.  From  December 
to  June,  western  part  and  eastern  part  of  the  basin  stay  nearly  constant  with  a  salinity 
value  of  almost  37  psu  and  almost  39.2  psu,  respectively.  With  the  beginning  of  July, 
eastern  part  of  the  basin  becomes  much  saline,  while  the  western  part  stays  nearly  stable. 
Maximum  salinity  is  observed  in  October  with  a  salinity  value  of  almost  39.8  psu  at  the 
eastern  part  of  the  basin. 

Figure  33.  Climatological  Monthly  Mean  Salinity  (psu)  at  the  Surface 
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Figure  34  shows  the  climatological  monthly  mean  salinity  at  150m  depth.  The 
winter  and  spring  seasons  show  a  similar  salinity  distribution  pattern.  Winter  and  spring 
seasons  are  much  saline  with  a  salinity  value  of  almost  39.2  psu  at  the  150m  depth 
comparing  to  the  winter  and  spring  seasons  at  the  surface  with  a  salinity  value  of  nearly 
38.7  psu.  Salinity  variation  shows  almost  a  stable  pattern  through  the  all  months  except 
the  Strait  of  Gibraltar  and  western  coasts  of  Africa  for  March  and  April. 


Figure  34.  Climatological  Monthly  Mean  Salinity  (psu)  at  150m 

Depth 
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Figure  35  shows  the  climatological  monthly  mean  salinity  at  600m  depth  with  no 
seasonal  variability  through  the  basin. 


Figure  35.  Climatological  Monthly  Mean  Salinity  (psu)  at  600m 

Depth 
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C.  ABSOLUTE  GEOSTROPHIC  VELOCITY  FIELD 

Many  studies  related  to  the  circulation  patterns  of  the  Eastern  Mediterranean  Sea 
have  been  explained  since  beginning  of  the  1900s.  Although  there  are  many  studies 
towards  the  Eastern  Mediterranean  Sea  circulation,  there  are  still  some  debates  about  the 
actual  circulation  patterns  of  the  area  (Alhammoud  et  ah,  2005).  The  structure  of  the 
actual  circulation  pattern  of  the  Levantine  Basin  in  the  Eastern  Mediterranean  is  way 
complex  and  it  is  needed  to  study  further  to  be  able  to  understand  the  actual  circulation 
patterns.  The  main  structures  are  basin-scale,  sub-basin  scale  and  mesoscale  structures. 
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Moreover,  there  are  permanent,  reeurrent  and  transient  eyelonic  and  anticyelonic  gyres 
and  eddies  in  the  interest  area  (Alhammoud  et  ah,  2005).  One  of  the  first  historical 
studies  was  carried  out  by  Nielsen  (1912).  Until  1960s,  the  knowledge  of  the  circulation 
pattern  of  the  Eastern  Mediterranean  Sea  was  limited.  Nielsen  explained  that  the 
circulation  was  a  basin-wide  cyclonic  gyre  and  the  strongest  current  was  close  to  the 
coast.  He  described  that  the  Coriolis  force  affected  the  currents  along  the  basin  coast 
(Alhammoud  et  ah,  2005).  A  new  general  circulation  schematic  was  provided  by 
Ovchinnikov  (1966).  Some  cyclonic  and  anticyelonic  sub-basin  scales  were  detected  by 
Ovchinnikov  (1966),  and  Lacombe  and  Tchernia  (1972)  in  the  1960s  and  1970s 
(Alhammoud  et  ah,  2005). 

In  the  Levantine  basin  (LB),  they  found  an  anticyelonic  gyre  at  the 
southeastern  corner  of  the  basin  and  a  cyclonic  gyre  between  Rhodes  and 
Cyprus,  which  might  be  related  to  the  Shikmona  anticyclone  (SMA)  and 
the  Rhodes  gyre  (RG),  respectively.  Ovchinnikov  (1966)  also  showed  that 
the  AW  flows  further  eastward  through  the  southern  part  of  the  Cretan 
passage  and  appears  to  reach  as  far  as  the  eastern  coast  of  the  LB.  Then 
this  AW  bifurcates  near  the  Egyptian  continental  shelf  One  branch 
continues  eastward  to  reach  the  eastern  edge  of  the  LB.  Then  it  turns 
northward  to  cross  the  Lattakia  Basin  up  to  the  northeast  of  Cyprus  and 
there  it  flows  westward  along  the  Asia  Minor  coast  and  between  the  Greek 
islands  to  exit  via  the  northern  part  of  the  Cretan  passage.  The  second 
branch  flows  offshore  to  join  its  twin  northwest  of  Cyprus  and  there  it 
merges  with  the  first  branch  described  above.  (Alhammoud  et  ah,  2005) 

After  the  Ovchinnikov’s  schematic  of  the  Mediterranean  Sea  circulation,  POEM 
group  made  surveys  in  the  Mediterranean  Sea  and  explained  a  different  pattern  from  the 
historical  circulation  schematics  of  the  Mediterranean  Sea. 

According  to  The  POEM  group  (1992),  among  these  changes  was  the 
presence  of  high  energy  and  variable  mesoscale  structures  consisting  of 
eddies  linked  by  meanders  and  jets.  The  POEM  group  (1992)  provided  a 
scheme  of  the  circulation  in  the  LB  characterized  by  two  systems  of  sub¬ 
basin  gyres,  an  anticyelonic  system  in  the  south  (Mersa-Matruh  and 
Shikmona  anticyclones;  MMA  and  SMA,  respectively)  and  a  cyclonic  one 
in  the  north  (RG).  The  two  systems  are  separated  by  a  strong  jet  called  the 
Central  Levantine  Basin  Current  (CLBC)  by  Ozsoy  et  al.  (1989)  or  the 
Mid-Mediterranean  Jet  (MMJ)  by  Robinson  et  al.  (1991).  This  jet,  which 
is  associated  with  the  AW  (characterized  by  its  salinity  minimum  values), 
enters  the  LB  via  the  Cretan  passage  and  meanders  eastward  between  the 
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RG  and  the  MMA.  Then  it  bifureates  south  of  Cyprus  into  two  branehes, 
one  merging  with  the  SMA  and  the  other  turning  northward  then  westward 
to  feed  the  Asia  Minor  Current  (AMC).  The  AMC  flows  westward  into  the 
Aegean  Sea  earrying  with  it  the  LSW  marked  by  its  high  salinity  and 
temperature.  (Alhammoud  et  ah,  2005) 

After  the  studies  of  the  POEM  group,  Mediterranean  Foreeast  System  (MFS) 
program  was  run  and  earried  observations  between  1998  and  2000.  The  result  of  this 
program  also  eonfirmed  that  sub-basin  antieyelonie  gyres  in  the  southern  part  of  the  basin 
were  persistent  (Alhammoud  et  ah,  2005).  Figure  36  shows  the  historieal  studies  about 
the  eirculation  patterns  of  the  Mediterranean  Sea.  Figure  36  (a)  is  from  historical  data  by 
Nielsen,  Figure  36  (b)  is  from  surface  geostrophic  current  by  Ovchinnikov,  Figure  36  (c) 
is  from  observations  between  1985  and  1995  by  POEM  group,  and  Figure  36  (d)  is  from 
XBT  observations  between  1998  and  2000  by  the  MFS  project.  Figure  37  shows  the 
circulation  features  schematic  of  the  Mediterranean  Sea  for  winter  and  summer  seasons. 
Figure  38  shows  the  circulation  features  of  the  Eastern  Mediterranean  Sea. 
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Figure  36.  Historical  Studies  of  the  Mediterranean  Sea  Surface 
Circulation  Pattern 


Source:  Alhammoud,  B.,  K.  Beranger,  L.  Mortier,  M.  Crepon,  and  1.  Dekeyser,  2005: 
Surface  circulation  of  the  Levantine  Basin:  comparison  of  model  results  with 
observations.  Prog.Oceanogr.,  66,  299-320. 
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Figure  37.  Circulation  Features  of  the  Mediterranean  Sea 


0  Llgurlan-Provencal  current 

1  Lion  Gyre 

2  Thyrrhenian  cyclonic  circulation  with  summer  weakening  and  eastern  anticyclone 

3  Algerian  Current  (3a)  and  eddies.  Atlantic-Ionian  Stream  (3e)  branching  of  the  Atlantic-lonian 
Stream  (3b.  3c.  3d)  and  Mid-Med.  Jet  (31) 

4  Rhodes  Gyre 

5  Western  Cretan  Gyre  12  Southern  Adriatic  gyre 

6  Western  Ionian  Gyre  13  Western  Adriatic  Co^ast  al  Current 

7  Anticylone  in  the  Gulf  of  Syrte 

8  Shikmona  and  Mersa-Matruh  gyres  system 

9  Cilician  and  Asia  Minor  Current 

10  lera-Petra  gyre 

1 1  Pelops  gyre 

Source:  Pinardi,  N.,  E.  Masetti,  2000:  Variability  of  the  large  scale  general  circulation  of 
the  Mediterranean  Sea  from  observations  and  modelling:  a  review.  Palaeogeogr., 
PalaeoclimatoL,  Palaeoecol.,  158,  153-173. 


AIS  =  Atlanoc-lonian  Stream  lA  =  Ionian  Anticyclones  PA  =  Pelops  Anticyclone 

MU  =  Mid  Ionian  Jet  MMJ  =  Mid-Meditefranean  Jet  CC  =  Cretan  Cyclone 

MAW  =  Modified  Atlantic  Water  ASW  =  Adriatic  Surface  Water  AMC  =  Asia  Minor  Current 

ISW  Ionian  Surface  Water  LSW  =  Levantine  Surface  Water 


Source:  Robinson,  A.  R.,  W.  G.  Leslie,  A.  Theocharis,  and  A.  Lascaratos,  2001: 
Mediterranean  sea  circulation.  Ocean  Currents:  A  Derivative  of  the  Encyclopedia  of 
Ocean  Sciences,  1689-1705,  doi:10.1006/rwos.2001.0376. 
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Figure  39  shows  the  upper  thermocline  cireulation  features  of  the  Eastern 
Mediterranean  Sea.  It  explains  the  various  observation  eruises  that  were  condueted  about 
the  Eastern  Mediterranean  Sea  cireulation  patterns  and  their  results.  Circulation  features 
have  some  characteristics  as  permanent,  transient,  and  recurrent.  Type  P  represents  the 
permanent  features,  type  R  represents  the  recurrent  features  and  type  T  represents  the 
transient  features.  Cyclonic  features  are  represented  as  C  and,  anticyclonic  features  are 
represented  as  AC. 


Eigure  39.  Upper  Thermocline  Circulation  Eeatures  of  the  Eastern 

Mediterranean 
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Source:  Source:  Robinson,  A.  R.,  W.  G.  Leslie,  A.  Theocharis,  and  A.  Lascaratos,  2001: 
Mediterranean  sea  circulation.  Ocean  Currents:  A  Derivative  of  the  Encyclopedia  of 
Ocean  Sciences,  1689-1705,  doi:10.1006/rwos. 2001. 0376. 


The  results  of  total  mean  absolute  geostrophic  velocity  and  climatological 
monthly  mean  absolute  geostrophic  velocities  at  different  depths  from  1960  to  2013  are 
seen  in  Eigures  40,  41,  42,  and  43,  respectively. 
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Figure  40.  Total  Mean  Absolute  Geostrophic  Veloeity  (cm/s)  at 
Various  Depths  Calculated  from  the  (T,  S)  Data  Using  the  P-Vector 

Method 
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Figure  41 .  Climatological  Monthly  Mean  Absolute  Geostrophic 
Velocity  (cm/s)  at  the  Surface  Calculated  from  the  (T,  S)  Data  Using  the 

P- Vector  Method 
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Figure  42.  Climatological  Monthly  Mean  Absolute  Geostrophic 
Velocity  (cm/s)  at  150m  Depth  Calculated  from  the  (T,  S)  Data  Using  the 

P- Vector  Method 


57 


Latitude  Latitude  Latitude  Latitude 


Figure  43.  Climatological  Monthly  Mean  Absolute  Geostrophic 
Velocity  (cm/s)  at  600m  Depth  Calculated  from  the  (T,  S)  Data  Using  the 

P- Vector  Method 
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The  results  of  the  identified  circulation  features  from  the  absolute  geostrophic 
velocities  of  the  Eastern  Mediterranean  Sea  are  listed  in  Table  7  for  0-1 50m  depth,  in 
Table  8  for  150-600m  depth,  and  in  Table  9  for  depths  below  600m. 
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Table  7.  Monthly  Mean  Cireulation  Features  in  the  Eastern  Mediterranean  Sea  at  0-150m  Depth 
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Table  8.  Monthly  Mean  Circulation  Features  in  the  Eastern  Mediterranean  Sea  at  150-600m  Depth 
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Table  9.  Monthly  Mean  Cireulation  Features  in  the  Eastern  Mediterranean  Sea  at  Depths  below  600m 
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It  is  observed  that  the  eireulation  features  weakens  as  depth  goes  down.  Some  of 
the  features  are  still  existent  in  some  months  at  the  150-600m  depth,  whereas  the 
majority  of  the  features  are  not  existent  at  the  depths  below  600m.  However,  there  might 
be  some  other  features  indueed  by  ehanging  atmospherie  faetors  and  thermohaline 
foreing. 

D,  HEAT  CONTENT 

1.  Surface  Layer  (0  -  150m) 

Figure  44  shows  the  monthly  variability  of  heat  eontent  at  the  surfaee  layer.  The 
heat  eontent  variability  is  nearly  zero  with  exeeption  of  Gulf  of  Gabes  in  Deeember. 
From  Deeember  to  April,  variability  is  negative  with  values  of  between  -0.4x10  J/m  and 
-1.4x10  J/m  .  June  shows  transition  eharaoteristies  between  negative  and  positive  heat 
eontent  variabilities,  while  May  shows  slightly  negative  heat  eontent  variability.  From 
beginning  of  June  to  September,  heat  eontent  variability  shows  a  positive  pattern.  Heat 
content  variability  increases  from  0.2x10  J/m  to  1.6x10  J/m  between  June  and 
September.  Heat  content  at  the  surface  layer  explicitly  shows  seasonal  variability. 
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Figure  44.  Monthly  Variability  of  Heat  Content  (J/m  )  at  the  Surface 

Layer 
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2,  Intermediate  Layer  (150  -  600m) 

Figure  45  shows  the  monthly  variability  of  heat  content  at  the  intermediate  layer. 
At  a  quick  glance  to  heat  content  variability  for  all  months,  it  is  observed  that  heat 
content  variability  is  slightly  up  or  down  on  the  zero  contour  line,  which  means  that  the 
heat  content  difference  is  almost  zero  with  exceptions  in  some  specific  areas  such  as  at 
the  comer  between  Turkey  and  Syria,  and  at  around  the  Gulf  of  Sidra  in  Libya.  There  is 
no  observed  significant  heat  content  variability  at  the  intermediate  layer. 
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Figure  45.  Monthly  Variability  of  Fleat  Content  (J/m  )  at  the 

Intermediate  Layer 
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3,  Deep  Layer  (600m  -  Bottom) 

Figure  46  shows  the  monthly  variability  of  heat  content  at  the  deep  layer.  The 
entire  basin  almost  shows  a  negative  heat  content  variability  pattern  with  some 
fluctuations  in  February,  June  and  July. 

When  the  three  different  layers  are  compared  at  once,  it  is  seen  that  the  surface 
layer  is  the  most  sensitive  one  to  the  heat  content  variability  that  is  affected  by  the 
external  mechanisms  such  as  winds,  solar  heating,  and  winter  cooling.  As  the  depth 
increases  through  the  intermediate  and  deep  layers,  monthly  heat  content  and  annual  heat 
content  nearly  balances  each  other.  Thus,  very  slight  variation  in  the  heat  content  is 
observed  at  the  both  intermediate  and  deep  layers.  It  is  observed  that  there  is  almost  no 
seasonal  heat  content  variability  at  the  deep  layer. 
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Figure  46.  Monthly  Variability  of  Heat  Content  (J/m  )  at  the  Deep 

Layer 
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E.  FRESHWATER  CONTENT 

1,  Surface  Layer  (0  -  150m) 

Freshwater  content  is  obtained  by  taking  the  difference  of  the  present  and 
reference  salinities  and  divided  by  the  reference  salinity,  which  is  40  psu.  Then,  it  is 
integrated  through  the  layer  depth.  After  the  integration,  freshwater  content  is  obtained  in 
meter  unit.  Figure  47  shows  the  monthly  variability  of  freshwater  content  at  the  surface 
layer. 
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Figure  47.  Monthly  Variability  of  Freshwater  Content  (m)  at  the 

Surface  Layer 
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From  January  to  June,  the  variability  of  freshwater  content  is  in  slightly  positive 
values  almost  through  the  basin  and  in  much  positive  values  in  the  entrance  of  Adriatic 
Sea,  and  at  the  Strait  of  Sicily.  From  July  to  December,  it  seems  that  the  variability  of 
freshwater  content  shifted  reversely  from  positive  to  negative  values  roughly  for  the  same 
locations. 

The  maximum  freshwater  content  variability  is  observed  in  March  and  April  at 
the  Strait  of  Sicily,  whereas  the  minimum  freshwater  content  variability  is  observed  in 
October  at  the  Gulf  of  Gabes,  which  is  located  at  the  coast  of  Tunisia  and  in  December  at 
the  Strait  of  Sicily. 
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2,  Intermediate  Layer  (150  -  600m) 

Figure  48  shows  the  monthly  variability  of  freshwater  content  at  the  intermediate 
layer.  It  is  observed  that  nearly  entire  basin  shows  a  negative  freshwater  content 
variability  for  all  months  except  August  with  a  very  minor  location  around  the  Tripoli. 
The  general  pattern  expresses  that  the  monthly  mean  freshwater  content  falls  behind  the 
annual  mean  freshwater  content. 

Figure  48.  Monthly  Variability  of  Freshwater  Content  (m)  at  the 

Intermediate  Layer 
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3,  Deep  Layer  (600m  -  Bottom) 

Figure  49  shows  the  monthly  variability  of  freshwater  content  at  the  deep  layer.  It 
is  observed  that  there  is  not  a  stable  pattern  of  freshwater  content  monthly  variability 
through  the  entire  basin.  July  draws  attention  about  high  positive  freshwater  content 
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variability  in  the  Ionian  Sea,  whereas  May  and  June  shows  a  high  negative  freshwater 
eontent  variability  in  the  same  area. 

Figure  49.  Monthly  Variability  of  Freshwater  Content  (m)  at  the  Deep 

Layer 

Monthly  Variability  of  Freshwater  Content  (600m-Bottom) 
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VI.  INTERANNUAL  VARIABILITY 


Empirical  Orthogonal  Function  (EOF)  method  is  used  in  this  ehapter  to  be  able  to 
analyze  the  inter-annual  variability  of  heat  content  and  freshwater  eontent  in  three  layers. 
Eaeh  mode  represents  a  percentage  of  variance.  In  this  study,  we  focus  on  the  first  three 
modes,  which  totally  cover  almost  45%  of  the  variances.  Inter-annual  variability  focuses 
on  anomalies  that  oecur  from  one  year  to  another.  So,  it  provides  information  about 
yearly  ehanges  for  a  long  time-series.  In  this  study,  inter-annual  variability  is  searehed  for 
54  years,  from  1960  to  2013. 

A.  FIRST  THREE  EOF  MODES  OF  HEAT  CONTENT  ANOMALIES  FOR 

THREE  LAYERS 

EOF  mode-1  covers  nearly  28%  of  the  time  ehanges  of  variance  for  the  Eastern 
Mediterranean  Sea.  Figure  50  shows  the  EOF  mode-1  and  assoeiated  Principal 
Component-1  (PC-1)  for  three  layers.  EOF  mode-1  represents  the  spatial  variation,  while 
the  PC-1  represents  the  temporal  variation  and  related  amplitudes.  Aceording  to  the 
temporal  amplitude,  it  ean  be  interpreted  that  in  whieh  area’s  heat  eontent  strengthens, 
while  the  other  area’s  heat  content  weakens.  At  the  surfaee  layer,  the  maximum  anomaly 
of  heat  content  is  observed  around  the  Crete.  The  PC  amplitude  maximum  is  observed 
nearly  in  1993  with  a  negative  sign,  whereas  the  PC  amplitude  maximum  is  seen  nearly 
in  2012  with  a  positive  sign.  Heat  content  anomaly  in  a  point  is  determined  by 
multiplying  the  EOF  values  and  PC  amplitudes.  For  instance,  heat  eontent  around  the 
Crete  in  1993  decreases  with  a  value  of  (0.03)  X  (-20  x  lO^J/m^)  =  -0.6  x  lO^J/m^. 
However,  same  spot’s  heat  content  in  2012  increases  with  a  value  of  (0.03)  X  (+20  x 
10  J/m  )  =  +0.6  X  10  J/m  .  The  red  eontour  line  shows  the  zero  magnitude,  where  there  is 
almost  no  heat  eontent  anomaly  for  the  areas  around  that  line.  For  temporal  pattern,  heat 
content  anomaly  is  positive  from  1960  to  1983  with  some  exceptions  such  as  1966,  1968, 
and  1977.  Positive  heat  content  anomaly  means  that  heat  eontent  inereases  within  that 
period.  From  1980  to  2000,  heat  content  anomaly  is  negative,  which  means  that  heat 
eontent  decreases.  From  2000  to  2013,  heat  content  anomaly  is  positive  again,  whieh 
means  that  heat  eontent  inereases. 
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At  the  intermediate  layer,  the  peak  heat  content  anomaly  occurs  around  the  Crete 
with  extended  coverage.  The  area,  where  peak  heat  content  anomaly  occurs,  runs  from 
south-west  coast  of  Greece  to  east  of  Crete.  In  this  area,  the  darkest  red  contour  expresses 
that  the  magnitude  of  the  heat  content  anomaly  is  larger  than  the  same  area  at  the  surface 
layer.  At  the  deep  layer,  it  is  observed  that  there  is  almost  no  heat  content  anomaly 
through  the  years. 


Figure  50.  EOF  Mode-1  and  PC-1  of  Surface  Layer  (0-150m), 
Intermediate  Layer  (I50-600m)  and  Deep  Layer  (600-bottom)  Heat 

Content  (J/m^) 


EOF  of  Heat  Content  (Mode:  1)  27.8% 


0.045 

0.04 

0.035 

0.03 

0  025 

0.02 

0.015 

0.01 

0.005 

-0.005 

-0,01 

-0015 

-0.02 

-0025 

-003 

-0.035 

-0.04 

-0.045 


Longitude 


Year 


70 


Figure  51  and  Figure  52  show  the  EOF  mode-2  and  EOF  mode-3,  respectively. 
The  idea  about  defining  the  heat  content  anomaly  behind  these  modes  is  same  with  the 
EOF  mode-1.  These  modes  represent  the  less  coverage  of  the  time  change  of  variances. 
The  surface  layer  in  the  EOF  mode -2  shows  almost  an  opposite  characteristic  with  the 
surface  layer  in  the  EOF  mode-1.  The  EOF  mode-2  heat  content  anomaly  at  the  surface 
layer  shows  negative  trend,  whereas  the  EOF  mode-1  heat  content  anomaly  at  the  surface 
layer  shows  positive  trend. 


Figure  5 1 .  EOF  Mode -2  and  PC-2  of  Surface  Layer  (0-150m), 
Intermediate  Layer  (150-600m)  and  Deep  Layer  (600-bottom)  Heat 

Content  (J/m^) 
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Figure  52.  EOF  Mode -3  and  PC -3  of  Surface  Layer  (0-150m), 
Intermediate  Layer  (150-600m)  and  Deep  Layer  (600-bottom)  Heat 

Content  (J/m^) 
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B,  FIRST  THREE  EOF  MODES  OF  FRESHWATER  CONTENT 

ANOMALIES  FOR  THREE  LAYERS 

EOF  mode-1  covers  nearly  25%  of  the  time  changes  of  variance  for  the  Eastern 
Mediterranean  Sea.  Figure  53  shows  the  EOF  mode-1  and  Principal  Component-1  (PC-1) 
for  three  layers.  EOF  mode-1  represents  the  spatial  variation,  while  the  PC-1  represents 
the  temporal  variation  and  related  amplitudes. 

At  the  surface  layer,  the  maximum  anomaly  of  freshwater  content  is  observed  in 
the  Ionian  Sea,  between  south-west  coast  of  Greece  and  Malta.  At  a  quick  glance,  the 
temporal  amplitude  is  observed  as  maximum  nearly  in  2013  with  a  negative  sign, 
whereas  the  temporal  amplitude  is  observed  as  maximum  nearly  in  1972  with  a  positive 
sign.  Freshwater  content  anomaly  in  a  point  is  determined  by  multiplying  the  EOF  values 
and  PC  amplitudes.  For  instance,  freshwater  content  around  the  Sparta  in  2013  decreases 
with  a  value  of  (0.015)  X  (-30  m)  =  -0.45m.  However,  same  spot’s  freshwater  content  in 
1972  increases  with  a  value  of  (0.015)  X  (-l-60m)  =  -l-0.9m.  The  red  contour  line  shows 
the  zero  magnitude,  where  there  is  almost  no  freshwater  content  anomaly  for  the  areas 
around  that  line.  For  temporal  pattern  until  2005,  freshwater  content  fluctuates  between 
positive  and  negative  values  with  some  positive  exceptions  such  as  1972,  1974,  and  1992. 
Positive  freshwater  content  anomaly  means  that  freshwater  content  increases  within  that 
period.  From  2005  to  2013,  freshwater  content  anomaly  is  negative,  which  means  that 
freshwater  content  decreases  within  that  period.  At  the  intermediate  layer,  the  peak 
freshwater  content  anomaly  develops  in  the  middle  of  the  Ionian  Sea.  West  of  the  Cyprus 
is  another  significant  spot  that  freshwater  content  anomaly  is  important.  In  this  area,  the 
darkest  red  contour  expresses  that  the  magnitude  of  the  freshwater  content  is  larger  than 
the  same  area  for  the  surface  layer.  At  the  deep  layer,  it  is  observed  that  magnitude  of  the 
freshwater  content  increases  in  the  core  of  Ionian  Sea  and  in  the  center  of  the  West 
Cyprus’  contour.  Freshwater  content  anomaly  at  deep  layer  is  larger  than  both  surface 
layer  and  intermediate  layers.  Figure  54  and  Figure  55  show  the  EOF  mode -2  and  EOF 
mode-3,  respectively.  The  idea  about  defining  the  freshwater  content  anomaly  behind 
these  modes  is  same  with  the  EOF  mode-F  These  modes  represent  the  less  coverage  of 
the  time  change  of  variances. 
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In  the  Levantine  basin,  freshwater  content  anomaly  at  the  deep  layer  is  much 
positive  than  both  at  the  intermediate  and  at  the  surface  layers’  freshwater  content 
anomalies  for  EOF  mode-1,  whereas  the  structure  is  reverse  for  EOF  mode -2. 


Figure  53.  EOF  Mode-1  and  PC-1  of  Surface  Layer  (0-150m), 
Intermediate  Layer  (150-600m)  and  Deep  Layer  (600-bottom)  Freshwater 

Content  (m) 
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PC  (m)  Latitude  {600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


Figure  54.  EOF  Mode -2  and  PC -2  of  Surface  Layer  (0-150m), 
Intermediate  Layer  (150-600m)  and  Deep  Layer  (600-bottom)  Lreshwater 

Content  (m) 
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PC  (m)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


Figure  55.  EOF  Mode -3  and  PC -3  of  Surface  Layer  (0-150m), 
Intermediate  Layer  (150-600m)  and  Deep  Layer  (600-bottom)  Lreshwater 

Content  (m) 
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VII.  CONCLUSIONS  AND  FUTURE  RECOMMENDATIONS 


A,  CONCLUSIONS 

The  data,  which  were  used  in  this  thesis,  were  obtained  from  World  Ocean 
Database  2013  (WOD13).  The  optimal  spectral  decomposition  (OSD)  method  was  used 
to  establish  the  synoptic  monthly  varying  three-dimensional  gridded  temperature  and 
salinity  data  with  0.25°x0.25°  horizontal  resolution  and  57  vertical  depth  levels.  The 
analysis  included  164,906  temperature  profiles  and  53,606  salinity  profiles.  After  the 
establishment  of  gridded  data,  the  seasonal  and  inter-annual  variability  of  thermohaline 
structure  and  circulation  of  the  Eastern  Mediterranean  Sea  were  investigated  from 
January  1960  to  December  2013. 

Sea  surface  temperature  increases  from  west  to  east  and  decreases  from  south  to 
north.  Sea  surface  shows  high  seasonal  temperature  variability,  whereas  the  depths  below 
surface  show  slight  seasonal  variation,  relatively.  Sea  surface  salinity  increases  from  west 
to  east.  There  is  almost  no  seasonal  salinity  variability  in  winter  and  spring  seasons,  while 
the  summer  and  fall  seasons  show  slight  seasonal  variation  at  the  surface.  Below  surface, 
there  is  no  seasonal  variation.  Sea  surface  absolute  geostrophic  velocity  shows  high 
seasonal  variability  in  circulation  features.  Below  the  surface,  absolute  geostrophic 
velocity  weakens  and  features  change  their  structures  along  with  their  magnitudes. 

EOF  analyses  were  studied  for  inter-annual  variability  of  heat  content  and 
freshwater  content  anomalies.  For  heat  content,  the  intermediate  layer  heat  content 
anomaly  is  larger  than  the  surface  layer.  However,  inter-annual  variability  of  heat  content 
anomaly  at  the  deep  layer  is  very  weak.  For  freshwater  content,  the  deep  layer  and 
intermediate  layer  show  similar  freshwater  content  anomalies  and  both  of  them  are  larger 
than  the  surface  layer.  It  is  observed  that  heat  content  anomalies  are  much  evident  and 
noteworthy  than  freshwater  content  anomalies. 
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B,  FUTURE  RECOMMENDATIONS 

The  Western  Mediterranean  Sea  is  the  eounter-part  of  the  Eastern  Mediterranean 
Sea.  Atlantic  Water  (AW)  provides  freshwater  flux  and  heat  flux  to  the  Eastern 
Mediterranean  Sea  through  the  Strait  of  Gibraltar.  The  circulation  patterns,  salinity,  and 
temperature  variations  in  the  Eastern  Mediterranean  Sea  are  widely  affected  by  the 
Western  Mediterranean  Sea.  Thus,  the  Western  Mediterranean  Sea  should  be  studied  at 
the  same  time  with  the  Eastern  Mediterranean  Sea  to  be  able  to  infer  much  clear  and 
detailed  result. 

In  this  thesis,  observational  data  was  used.  The  research  of  ocean  variability 
depends  on  both  studies  with  observational  data  and  modeling  studies.  Sometimes,  the 
observational  datasets  may  be  insufficient  or  separated  temporally  and  spatially.  Thus,  it 
may  not  allow  us  to  make  detailed  research  about  the  mechanisms  that  drive  the 
variability  of  the  Eastern  Mediterranean  Sea.  So,  a  modeling  study  along  with  the 
observational  data  studies  may  be  useful  in  order  to  have  a  better  understanding. 

The  Adriatic  Sea  and  Aegean  Sea  are  the  other  components  of  the  Eastern 
Mediterranean  Sea.  Eor  better  results,  those  seas  should  be  included  in  the  future 
researches. 

Ereshwater  content  anomaly  shows  a  negative  pattern,  while  the  heat  content 
anomaly  shows  a  positive  pattern  since  2005.  This  might  be  due  to  the  global  warming 
and  climate  change.  Thus,  a  further  study  may  be  useful  to  be  able  to  interpret  this 
relation. 
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APPENDIX  A.  MONTHLY  MEAN  TEMPERATURE  AT 
DIFFERENT  DEPTHS 
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Monthly  Mean  Temperature  at  Depth:  20m 
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Monthly  Mean  Temperature  at  Depth:  50m 
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Monthly  Mean  Temperature  at  Depth:  100m 
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Monthly  Mean  Temperature  at  Depth:  150m 
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Monthly  Mean  Temperature  at  Depth:  200m 
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Monthly  Mean  Temperature  at  Depth:  250m 
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Monthly  Mean  Temperature  at  Depth:  300m 
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Monthly  Mean  Temperature  at  Depth:  400m 
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Monthly  Mean  Temperature  at  Depth:  600m 
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APPENDIX  B.  MONTHLY  MEAN  SALINITY  AT  DIFFERENT 

DEPTHS 
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Monthly  Mean  Salinity  at  Depth:  10m 
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Monthly  Mean  Salinity  at  Depth:  20m 
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Monthly  Mean  Salinity  at  Depth:  100m 
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Monthly  Mean  Salinity  at  Depth:  150m 
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Monthly  Mean  Salinity  at  Depth:  200m 
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Monthly  Mean  Salinity  at  Depth:  250m 
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Monthly  Mean  Salinity  at  Depth:  300m 
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Monthly  Mean  Salinity  at  Depth:  400m 
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Monthly  Mean  Salinity  at  Depth:  600m 
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Monthly  Mean  Salinity  at  Depth:  1000m 
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APPENDIX  C.  MONTHLY  MEAN  ABSOLUTE  GEOSTROPHIC 
VELOCITY  AT  DIFFERENT  DEPTHS 
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Monthly  Mean  at  Depth:  20m 
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Monthly  Mean  at  Depth:  50m 


10  15  20  25  30  35  10  15  20  25  30  35  10  15  20  25  30  35 

Longitude  Longitude  Longitude 


104 


Latitude  Latitude  Latitude  Latitude  Latitude  Latitude  Latitude  Latitude 


Monthly  Mean  at  Depth:  100m 
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Monthly  Mean  at  Depth:  300m 
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Monthly  Mean  at  Depth:  600m 
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PC  (10®  J/m^)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


APPENDIX  D.  FIRST  FIVE  EOF  MODES  OF  HEAT  CONTENT 


EOF  of  Heat  Content  (Mode:  1)  27.8% 
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PC  (10®  J/m^)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


EOF  of  Heat  Content  (Mode:  2)  11% 


40 

38 

36 

34 

32 

40 

38 

36 

34 

32 

40 

38 

36 

34 

32 


0.035 

0.03 

0.025 

0.02 

0.015 

0.01 

0.005 

-0.005 

-0.01 

-0.015 

-0.02 

-0.025 

-0.03 

-0.035 


10  15  20  25  30  35 


Longitude 


no 
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EOF  of  Heat  Content  (Mode;  3)  6.02% 
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PC(10®J/m^)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


EOF  of  Heat  Content  (Mode:  4)  4.65% 
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PC(10®J/m^)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


EOF  of  Heat  Content  (Mode;  5)  4.27% 
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PC  (m)  Latitude  (600ni-Bottom)  Latitude  (150m-600m)  Latitude  (0-1 50m) 


APPENDIX  E.  FIRST  FIVE  EOF  MODES  OF  FRESHWATER 

CONTENT 


EOF  of  Freshwater  Content  (Mode:  1)  25.2% 
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PC  (m)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


EOF  of  Freshwater  Content  (Mode:  2)  11.9% 
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PC  (m)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


EOF  of  Freshwater  Content  (Mode:  3)  8.47% 
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PC  (m)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


EOF  of  Freshwater  Content  (Mode;  4)  6.3% 
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PC  (m)  Latitude  (600m-Bottom)  Latitude  (150m-600m)  Latitude  (0-150m) 


EOF  of  Freshwater  Content  (Mode:  5)  5.65% 
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